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I. Introduction 

Direct evidence for the top quark was recently presented by the CDF (Collider 
Detector at Fermilab) Collaboration This is considered to be a remarkable suc- 
cess for the Standard Model (SM), since the value of the top mass determined by 
them, 174 ± 10tj2 G^^V, coincides with the indirect determination from the available 
precise data of electroweak experiments. Due to the determination of the top quark 
mass, it is certain that a number of new and interesting issues in top quark pair 
production and decay will arise. At the Next Linear Collider (NLC), operating in 
photon-photon collision mode at a design energy of 500 ~ 2000 GeV with a lumi- 
nosity of the order lO^^cm^^s^^ 0, a large number of heavy quark pairs and other 
new particles can be produced with an agreeable production rate . Those events 
would be much cleaner than those produced at pp and pp colliders. It has been 
also found that the ti production rate in 77 collisions from laser back-scattering is 
much larger than that from the direct e~^e~ —>■ ti production both with and without 
considering the threshold QCD effect, for ~ 100 — 20{)GeV at center-of-mass 
energies of the electron-positron system around 1 TeV.Q. Thus the process 77 — > ti 
has a large potential for producing and studying heavy quarks directly. The next- 
to-leading order QCD corrections in the SM for this process both for polarized and 
unpolarized particle collisions have been discussed in detail in Ref.0. There it was 
shown that these QCD corrections are significant, about 10%. However, the process 
also seems to be sensitive in searching for virtual effects from new particles, which 
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would imply physics beyond the Standard Model. Among various models which 
introduce new physics, the Minimal Supersymmetric Standard Model (MSSM)0] is 
the most attractive one at present, since it is the simplest case of Supersymmet- 
ric Models (SUSY). Therefore it is important to study the physical effects induced 
by the virtual supersymmetric particles for top quark pair production via photon 
fusion. If we assume that one of the stop quarks are much lighter than the other 
supersymmetric quarks due to strong mixing, the process 77 — > ti would be a very 
interesting reaction for testing supersymmetric theory in that the corrections due 
to virtual t would provide evidence for the existence of the quark scalar partners. 
Recently H. Wang el al. have calculated the Supersymmetric QCD corrections for 
this process via unpolarized photon collisions in some cases in Ref. . 

With the advent of new collider techniques, back-scattered laser photons can 
be linearly polarized with high luminosity and rather high efficiency These 
techniques will undoubtly be used to improve our knowledge of the top quark pa- 
rameters, particularly when one combines the data from top pair production via 
photon-photon fusion with other data from e'^e~ and pp collisions. In addition to 
that general advantage, the process of top quark pair produced by polarized photon 
collisions may provide us with a better understanding of CP-violation phenomena 
0. Since the mixing of top quarks with other generations is very small, CP- violation 
in this process is negligibly small within the Kobayashi-Maskawa framework of the 
Standard Model. However, in the MSSM, even without generation mixing, there 
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are more possibilities to introduce complex couplings than in the SM. Here in this 
paper, we limit ourselves only to the combination of the CP-violating phases 0^ and 
(f)g, which emerge in the stop quark mixing matrix and in the Majorana mass term 
of the gluino |Ty]. These non-zero complex phases, which occur in the Lagrangian, 
cannot be rotated away by a suitable redefinition of the fields. They will lead to CP 
violation within a single generation. Therefore, the observation of CP-violation in 
the process of top quark pair production via polarized 77 fusion may help to pro- 
mote our understanding of the features of the combination of CP-violating phases 
(pl and (j)g. 

The paper is organized as follows: in Sec. II, the renormalization of the top 
quark self-energy in the frame of the MSSM is described, and the explicit analytical 
forms of the Lorentz invariant matrix elements including the next-to-leading order 
of supersymmetric QCD corrections are presented. In Sec. Ill, the numerical results 
and discussions are depicted. Finally, the conclusions are given. In the appendix, 
the form factors used for the cross section calculations are listed. 

II. Calculations 

In this work, we denote the reaction discussed as 



7(P3, Ai)7(P4, A2) t{p2)i{vi) (2.1) 
where Ai,2 = ±1, P2 and pi represent the momenta of the outgoing top quark 
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and its anti-particle, ps and p4 are for the momenta of the two incoming photons 
with hehcities Ai and A2 respectively. 

The Feynman diagrams for the process 77 —>■ tt are shown in Fig.l, and the 
relevant Feynman rules can be found in [^]. In the calculation dimensional regu- 



larization and the on-mass-shell (OMS) renormalization scheme |[TT|] are adopted to 
eliminate all the ultraviolet divergences. The interaction Lagrangian for gluino-stop- 
top coupling including the CP-violating phases is given by 

+ {sin9i ii - cosOi t2)e"^PL] Qa + h.c. (2.2) 

= -V^QsT'^ t [{Ar h + Br t2)PR + {Al h - Bl i2)PL] 9a + h.C. 

where Qs is the strong coupling constant, T" are SU(?>)c generators and Prl = 
i(l ± 75). 6^ is the stop mixing angle which transforms the stop mass eigenstates 
tn (n=l,2) to the weak eigenstates and Ir. (j) = (pi — (pg is the combination of the 
two phase angles. In the Equation (2.2) we denote 

AR = coseie-'t Al = sin9ie''t>, 

Br = sinOie-'^, Bl = cosOie^'^ ^ ' 

With the introduction of the CP-violating phase in the MSSM, the renor- 

malized one-particle irreducible two-point function for top quarks containing the 



contributions from SUSY QCD should be written as |T2 



f(p) 



i{i> -mt) + i i)PL^^{p^) + ^PrJ:^{p^) + pj:^ (p^) + PrS* ^p 



(2.4) 
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It should be mentioned here that in the above equation the upper conjugation 
symbol | acts only on the CP-violating phase 0. The corresponding unrenormalized 
top quark self energies read 



S^(P') = -^gl(An,BLB,[p\mg,mi^]+ALBRBr\p\mg,miS). (2.5) 



^''(P') = iALBRB,[p\mg,mi^] + A^i?z.i?ib^ m^, m,-j) , (2.6) 



S^(p') = f^glm-,ALBL(B^[p\mg,mi^]-B^[p\mg,miS), (2.7) 



S^*(P') = ^glmgARBR[B^\p\mg,mi^]-B,[p\mg,miS)- (2.8) 

rrit and mi^ {n = 1,2) involved above are the masses of the top quark and stop 
mass eigenstates. In the SU(3) group, Cp = {N^ — l)/{2Nc). Imposing the on- 
shell renormalization conditions given in Ref. PT|, one can obtain equations for the 
renormalized self-energy functions: 



mtRe±^{mf) + Ret^ ^{ml) = 
mtRet^{ml) + Ret^{ml) = 
mtRet^iml) + Ret^iml) = 



mtRet^{ml) + Ret^ ^ml) = 

i?eE^(m,2)+i?eE^K)+2mi^i?e(mtE^(/)+miE^(p^)+E^(p^)+E^ = 

(2.9) 

where Re only takes the real part of the loop integral functions appearing in the 
self energies. The relevant self-energy and vertex counterterms are 



(2.10) 



tSE = zICl^Pl + Cr^Pr-CsPl-C+Pr], 

iSAf^ = -ieQ,^^^[C-PL + C+PR]. 
where Qt is the charge of the top quark and 

Cl = C- = USZ'^ + SZ'^^), 

Cr = C+ = UdZ^ + 5Z^^), 

C- = !f(SZ^ + SZ'^^) + Sm, ^ ' 

C+ = iflsz^ + 5Z^'') + Sm. 

We choose the field renormalization constant SZ^ to be real. Using Eq. (2.9)- 
(2.11), we obtain the renormalization constants as 



Sm = I (mti?eE^(m2) + mti?eE^(m2) + i?eE^(m?) + i?eE^ t(m2)) , (2.12) 



SZ^ = -E^(m2)-^ 

\ t / nit 



Rell^HmD-Rell^imV, 



mtReE^{p^) + mtReT^^if ) + Re^^{p^) + i?eE^ ^{p^) 
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(2.13) 



5Z^ = -E^(m2) -mt-^ [mti?eE^(p2) + mtReE^{p'^) + i?eE^(p2) + i?eE^ ^{p^) 

(2.14) 

Including all next-to-leading order supersymmetric QCD corrections, the renor- 
malized amplitude for tt pair production in 77 collisions is shown as 



\p^=m? 1 



M,e„(A2,Ai) = Mo(A2,Ai)+5Mi-'-f(A2,Ai) 

= Mo(A2, Ai) + 5M^^^/(A2, Ai) + 5M^^'-*^^(A2, Ai)+ (2.15) 
(5M''°^(A2,Ai) + (5M*'-(A2,Ai) 

where Mq is the amphtude at the tree level, SM^^^f, ^M''^''*^^, SM'""' and 

represent the renormalized amplitudes with the next-to-leading order supersymmet- 

ric QCD corrections arising from the self-energy, vertex, box and triangle diagrams, 

respectively. Their explicit forms are given by 



Mo(A2, Ai) = 



- ij^e,M{P4, A2)e^(p3, Ai)u(p2) l'' ih -i>A + rnt)Y viPi) 

(2.16) 



5M*(A2, Ai) = i-^ ^e^(P4, A2)e^(p3, Ai) «iv(P2) 

"^t ) N=L,R 
+ iP3 ^P4,IJ'^ i^, i u) 



(2.17) 



(5M^(A2,Ai) = -i f _ *2 ^m(P4, A2)e^(p3, Ai) ^ un{p2) 

^ '^t N=L,R 

iftS + /2,S K 7^ + fS p'2 Y 

+ Pi P2 + C ^^47^7-^ + i>^\l' + n% MY) v{p,) (2-18) 
+ (P3 ^Pa,I^^ i^, i u) 



5M^{\2,\i)^i ^e,.(p4, Ai)e,.(p3, A2) X! ^n{P2) 

N=L,R 
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+ /SIp^P^ + f'SU^l'r + /SUt'^T'^ + 04^5*^^ (2.19) 

+ (P3 ^ P4, 1^ ^ I'J ^ u) 

SM'^{X2, Ai) = A2)6.(P3, Ai) ^ /;^«iv(P2)^;(pi) (2.20) 

^ N=L,R 

Here £ = - p2)^ = (pi -Pif- The explicit form factors /^SJ> /w^ 

are presented in the appendix. The analytical deduction for the renormalized ampli- 
tudes in Eq.(2.15) clearly states the complete cancellation of ultraviolet divergences, 
which is required to ensure the correctness of our calculation. 

The cross section of the process for the polarized photons is given by 



where Nc is the number of colors, — (m| — \s) ± \s(3t , (5t — \J^ — 'im.t/s. 
Note that the summing over the spins is performed only over the final quark pair: 



Y: \M^i;^^{s,i)\' = Y: |Mo"-"f + 2Re[Y: M^^'^^ ■ SMttp I ■ (2.22) 

spins spins \ spins 

For polarized massless vector particles we have 
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where p is an arbitrary light-like Lorentz vector. 
The cross section with unpolarized photons is 



a = - ^(Ai,A2). (2.24) 

^ Ai,A2=-l 

III. Numerical results and discussion. 

In the numerical evaluation, we take the top quark mass to be nit = 175 GeV , 
a = 1/137.036 and use the two-loop running coupling constant a^. For the MSSM 
parameters, we assume that the supersymmetric weak eigenstate partners ii and 
mix equivalently, namely the mixing angle 6^ = 45°. For the other two independent 
parameters in stop quark mass matrix, we use the masses of the stop mass eigen- 
states, and assume m^^ < m^^- Current experiments |[TB| constrain the light stop 



mass eigenstate to be heavier than 50 GeV. To describe CP violation in tt production 
in polarized 77 fusion, a CP- violating parameter ^cp is defined as below: 

a{++) + a{ ) 

We set the CP-violating phase to be = 45° when evaluating the cross section 
in the case of non-zero phase angle, so that the CP- violating effect can be maximized. 

In Fig. 2(a) (b) the corrections for the cross sections cr(±±) and o"(±=l=) for dif- 
ferent initial photon linear polarizations are plotted as functions of the cm. energy 
V^, assuming the following mass values for supersymmetric particles involved: 
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m-g = 100 GeV, mi^ = 100 GeV, = 450 GeV (3.2) 

Unlike the corrections from the SM 0, all the contributions of the MSSM are 
negative in all the polarization cases. The relative correction 6, defined as {a — 
o"o) / (3"o, is generally in the order of 10~^ which is less than that of the Standard Model 
QCD corrections p. But it can also reach about 10% in some points. Fig. 2(a) and 
Fig. 2(b) show that the corrections will be suppressed in the case of introducing a 
non-zero CP-violating phase angle for both of the polarized photon fusion modes 
under the conditions of Eq.(3.2). The dependence of the CP- violation parameter 
^cp on \/§ is given in Fig.2(c). It shows that in the cm. energy range around the ti 
threshold, CP is strongly violated, where C,cp can be above 2.6%. The curve drops 
sharply with increasing V^. From the definition of the CP-violating parameter ^cp 
given in Eq.(3.1), one can understand that the feature of the curve in Fig.2(c) is 
simply the result from Fig.2(a,b). 

Since the present experimental results can not exclude the existence of very light 
gluinos and there has been renewed interest in this case recently, we also calculate the 
relative corrections to the cross sections with a very light gluino mass, as functions 
of when 

m-g = 5 GeV, mi^ = 50 GeV, mi^ = 250 GeV (3.3) 

The relative corrections S{±, ±) and 5(±, =i=) for different photon linear polariza- 
tions are plotted in Fig.3(a)and Fig.3(b) respectively. The curves go from positive 
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to negative values as \f§ increases. The relative corrections are positive around tt 
threshold, whereas they are always negative with the conditions Eq.(3.2) where the 
gluino mass is much heavier. Although in most energy regions the absolute correc- 
tions are small compared with the results given in Fig.2(a,b), the maximal value 
of relative correction is still about 3.5% around — 500 GeV which can be seen 
in Fig. 3(a). The high peaks shown in Fig.3(a) are the results from resonance effect 
when ~ 2mf2 = 500 GeV, which comes from the contributions of the triangle 
diagrams in Fig.l(i). The CP violation quantities in this case are plotted in Fig.3(c). 
All the values of ^cp have negative signs, and are about an order of magnitude lower 
than the values in Fig. 2(c), varying from —0.184% to —0.296%. 

As a check, we also calculate the corrections with extending cm. energy \/§ 
up to the region between 2 ~ 5 TeV and taking the same parameters of Eq.(3.2) 
and Eq.(3.3) respectively. The absolute corrections get smaller with increasing V^. 
Noteworthily the convergency of the correction with heavy gluino and squarks is 
much faster than that with light gluino and squarks. It proves that the decoupling 
really takes place at high energies, and the results of our cross section spectra at 
one-loop level satisfy the unitarity which is expected on theoretical grounds. 

According to Eq.(2.24), we can get the relationship between the unpolarized cross 
section and the gluino mass, which is shown in Fig.4(a) with 

rrii^ = 100 GeV, mi^ = 450 GeV and = 500 GeV (3.4) 
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From the figure, we can see that the corrections jump abruptly at the point rrig = 
75 GeV. This sudden change is due to the influence of a singularity in top quark wave 
function renormalization at the point rrit — rUg + m^^ . This singularity originates 
from the renormalization constants 5Z^ and 5Z^ of the top quark wave function 
shown in Eqs.(2.13) and (2.14). The effect of the singularity makes the absolute value 
of the correction get a large increment near the singularity point. However at the 
gluino mass range which is far away from the singularity mass region, the correction 
approaches to a small value quickly. This implies that there exists decoupling with 
heavy gluino in the SUSY QCD corrections. The corresponding prediction for the 
CP- violation parameter ^cp as a function of rrig is shown in Fig.4(b). There the 
singularity effect from the top quark wave function renormalization exists at the 
point rrig — 75 GeV . The inffuence of the singularity also can be seen from the 
fact that in Fig. 3(a) (b) the corrections are distinctly much smaller than those in 
Fig.2(a)(b). That is just because the parameters for Fig.3(a)(b) as shown in Eq.(3.3) 
are far away from the singularity point. 

The relations between the corrections of the unpolarized cross section and the 
two stop quark masses m^^ and m^^ are depicted in Fig.5(a) and Fig. 6(a) respectively. 
In Fig. 5(a), the large corrections near the point m^^ — 75 GeV are again due to the 
inffuence of the singularity. However, the corrections show the decoupling with the 
masses of stop quarks, when mf. (i = 1, 2) are far away from the singularity points 
mi. = rrit — nig. The decoupling effect can be seen more clearly in the heavy stop 
case. In Fig. 6 (a) the curves of corrections are rather ffat in the heavy t2 region, 
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and the small waves around m^^ = 250 GeV are merely due to the resonance effect 
when = 2m(2- The corresponding dependence of parameter ^cp on these two 
stop quark masses is given in Fig. 5(b), 6(b). 

Fig. 4 (a), Fig. 5 (a) and Fig. 6 (a) also show clearly that the absolute corrections 
for = 45° are apparently smaller than those for ^ = 0° in unpolarized photon 
fusion, except in the light gluino mass region in Fig.4(a) which is simply due to 
the singularity effect in top quark wave function renormalization. The fact that 
the corrections are always maximal for vanishing CP violation is simply due to the 
fact that with introducing a non-zero CP- violating phase angle in the gluino-stop-top 
coupling, the real part of this coupling strength will be reduced down quantitatively. 
Further analyses show that this reduction will lead a decrease in the absolute values 
of the real parts of form factors, whereas the imaginary parts of form factors do 
not contribute to the cross section. But around m^^ ~ V§/2 mass region shown in 
Fig. 5 (a), the situation is more complicated and this feature does not exist. 

IV. Conclusion 

In this work we have studied the one-loop supersymmetric QCD corrections of the 
process 7(A2)7(Ai) tt. 

From the results of numerical calculation, we can conclude that the singularity of 
the top quark wave function renormalization affects the corrections heavily. Large 
corrections can be expected near the singularity region which satisfies = rrig + 
mi. {i — 1,2). This effect is shown in Fig.4(a), 5(a) and 6(a), where the large 
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corrections are all due to the singularity effect. This singularity effect can be also 
used to explain why the corrections are large in Fig. 2 (a), whereas the corrections in 
Fig. 3 (a) are small where the parameters are far away from all the singulatriy points. 
However the corrections are insensitive to the masses of supersymmetric particles 
ti {i — 1,2) and g when their masses are heavy and far away from singularity points. 
This feature is just the manifestation of the decoupling effects. 

In any case, with the fusion of photons polarized with parallel spin directions, 
the corrections are always significant near the tt threshold. Furthermore, if the CP- 
violating phase really exists in the squark mixing matrix or in the Major ana mass 
term of the gluino predicted by the MSSM, the CP- violating parameter could be 
expected to be the order of 10~^ ~ 10~^. Therefore testing CP-violation in this 
process is one of the promising tasks for a future photon-photon collider. 

This work was supported in part by the National Natural Science Foundation 
of China and the National Committee of Science and Technology of China. Part of 
this work was done when two of the anthers. Ma Wen-Gan and Yu Zeng-Hui, visited 
the University Vienna under the exchange agreement (project number: IV.B.12). 



15 



Appendix 



The form factors //j^ can be expressed as 

s(t) 



/r,§ = ^(P2 ■ Pi) [ALBRfntBi [-P2 + Pi, m-g, m^-J + {L ^ R, 



4^{P2 ■ Pa)AlBl {mgBo[-p2 + P4, m-g, mjj + {m-g -m-g, m^^ m^J } 
+ '^{P2-P4)iCs -mtCn) 

(Al) 



f^R = ^ { C^^lBr P2-P4- ml)Bi [-P2 + Pi, m-g, m^J + {L ^ R, m^^ m^J } 

+ ^^{AlBl + AuBRjmt |mgSo[-p2 + Pi, m-g, mfj + (m^ -mg, mi^ m 
- 2mt{Cs+Ct) + 2mj{CL + CR)-A{p2-pi)CR 



where the definitions of the renormalization quantities Cf and Cr^l can be found 

in Eq.(2.11), and those of Ar^l and Br^l in Eq.(2.3). The replacement mg — > —mg 

is not performed on the argument of loop integral functions B, C and D. 

The form factors of triangle, vertex and box diagrams are given as follows 

jtr ^ nnc^ {{ArBlCii - ArBlCi2 + AlBrCi2) [-P2,Pi + P2, m-g, m^^ , m,-J 
+(L ^ R,mi^ mfj} 
- ^AlBl {m~gCQ[-p2,Pi +P2,mg, m^^, m^J 

{A5) 



16 



fiS = ^f^t {{AlBr - ArBl)C24[-P2,P4, rrig, mi^.mi^] + {L ^ R, ruf^ m^-j} 
+ mt{CR-CL) 

m 



& = -^2(pi • Pa) [AlBr{Ci2 + C23)[pi, -P3, m~g, mi^.mi^] + {L ^ R, mi^ m^j} 



fin = #^ {(4>liS«C24 - m2(Cn + C21) + 2 p2 • P4 ^^^^(Cis + C23))[-P2,P4, m^, m,-^ 
(L 4^ R,mi^ ^rni^)} 

+ 4^^ 

(^9) 



>lii-BLC23 - ALSijC23)[pi, -pa, mp, mj-^, mfj + (L ^ i?, mj^ mtj} 
- ^AlBl {mg(Co + Cii) [pi, -p3, TOg, rrii^ , m^J + (m^ -rrig, mi^ m^J } 

(AlO) 



/5,R = ^ {^L5i?C24[-P2,P4, rrig, mj^, rn^-J + (L ^ i?, m^^ ^ mj-j} 
+ 2Ci^ 



17 



Mt) 
J7,R 



'-^mt {{ArBlCu - ArBlC,2 + AlBrC,2 + ArBlC2i- 
ArBlC23 + AlBrC2z) [-P2,P4, rUg, mi^ , m^J + (L ^ R, mi^ J } 



-^AlBl {mg(Co + Cii)[-p2,P4, m~g, mi^.mi^^ + {m~g -rrig, ruf^ m^-j} 

(^13) 



{AU) 



Mt) 
Jl,R 



= ^mt {{ArBlD^u - ArBlD^^^ + ^^5^^313) 

[-P2 ,P4,P3,mg, mi^ , mi^ ,mi^] + {L ^ R, ) } 

- ^AlBl [mgD27[-p2,P4, P3, m-g, m^^ , m^^ , m^J + {m-g -m-g, mi^ 



m 



t2 



(A15) 



)} 



Mi) _ Mt) 



J2,R — J: 



l,R 



(A16) 



h,R = -^AlBr{D27 + Dzi2)[-P2,PA,P3,m-g,mi^,mi^,mi^] ^^-^^-^ 



/4,S = -^AlBrDsis [-P2,P4, P3, rrig, mi^ , mi^ , J 
+ {L ^ R,mi^ ^ rriij 



(AIS) 



/5S = ^AlBr{D3ii - L>3i2) [-P2,Pa, P3, m-g, mi^ , m^-^ , m^- J 
+ {L ^ R,mi^^ mi^) 



= ^AlBr{D27 + -D311 - -D313) [-P2, P4, P3, ^fl, ^ti , ^ti , ^tj 



{A20) 
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/tS = ^mt{{ARBLD23-Ai^BRD2s-ARBLD25-ARBi^Dsio+ 
ArBlDsq - AlBrD^q) [-P2, P4, Pi.rn-g, mi^ , m^^ , J 
+(L ^ R, mi^ mi^)] (A21) 
+ ^AlBl {mg(Di3 + D26) [-P2,P4,P3, rrig, mi^ , rrif^ , m^J 



m = {{ArBlDs, - ArBlDsw - ArBlD37 + AlBrDs7+ 

ArBlD^^ - AlBrD^^) [-P2, P4, P3, m~g, mi^ , m^^ , J 
+{L^R, mi^ mjj } (A22) 
- ^AlBl {mg{D25 - D2e) [-P2, P4, Ps, rUg, mi^ , mi^ , m^J 
+{m~g -m~g,mi^ m|-j} 



f9,R = ^rnt {{ArBlD^^ - ArBlD,^ + AlBrDi^ + ArBlD2i+ 

ArBlD23 - AlBrD23 + ArBlD2A - 2ArBlD25 + 

AlBrD25 - ArB^D2, + A^BnD^e - 2ArBlD^^q+ 
AlBrD^.o + ArBlDu + ArBlDsq - AlBrD^^) 
[-P2 ,P4,P3,mg, mi^ , mi^ , mi^ ] (^23) 

- ^AlBl {m-g{Do + + D^2 - D^s + D2A - D2q) 
[-P2,P4, P3, rrig, mi^ , mi^ , m^J 
+{m~g -rUg, mi^ m^j} 

^mt {{ArBlD24 - ArBlD21 + ArBlD25 - AlBrD25- 
ArBlD2q + AlBrD26 - ArBlD^i - 2ArBlD3w+ 

AlBrDsio + ArBlD^^ + 2ArBlD35 - AlBrD^^,- 
ArBlD^j + AlBrDs^ + ArBlDsq - AlBrDsq) 
[-P2,P4,P3,m~g, mi^ , mi^ , mj^ ] (^24) 
+ (L ^ R,mi^ m^j} 

[-P2, P4, P3, ^g, ^fi , , "^tj 



JlO,i? 



+ 
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Mi) _ 

Jll,R — 



-^AlBr{Dsi2 - D313) 



(A25) 



+ 



fKi) 
Jl2,R 



Mt) 
Jll,R 



{A26) 



Mi) 

Jl5,R 



Mi) 
JlQ,R 



where 



/S = ^AiSR(L>26 + ^38 - ^39 - i^23) 

/S = ^Az.5R(D37 + ^38-^39-i^310) 
[-P2, P4, P3, ^g, ^ti , , ^tj 

+ {L ^ R,mi^ ^ rriij 

= ^AlBr{Di3 - D,2 - D22 - D23 - D2A + i^25 + 
+ -D3IO ~ -D36 + -D38 ~ -D39) 

[-P2, P4, P3, mg, , , rnj-J 
+ {L ^ R,mi^^miJ 

= ^AlBr{D24-D22-D25 + D26 + D3^-Dss- 
D36 + L>37 + i^38 - ^39) 

[-P2 ,P4,P3,mg, mi^ , mi^ , rrif^ ] 



(P3 - Pi)^ = ^, s = (Pi + P2)^ 
2pi ■ P3 = 2p2 ■ P4 = ml - 
2p3 -Pi^s, 2pi ■p^^2p2-p3^ml - it 




2pi • P2 = s — 2mj , s + i + It = . 
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All the definitions of loop integral functions A, B, C and D used in our paper 



can be found in Ref.|Tj]. The numerical calculation of the vector and tensor loop 
integral functions can be traced back to the four scalar loop integrals Aq, Bq, Cq 



and Do in Ref. [0^5 
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Figure Captions 

Fig.l Feynman diagrams for contributions from the tree- level and next-to-leading 
order in Supersymmetric QCD for the 77 — > tt process: (a-b) tree level diagrams; 
(c-e) self-energy diagrams; (f-g) vertex diagrams; (h) box diagrams: (i) triangle 
diagrams. The dashed lines represent ^1,^2 for (c-i). 

Fig.2 With nig — m^^ — 100 GeV and m^^ — 450 GeV: (a) the cross section 
of the ti production process cr(±±) as a function of v^, solid line for tree-level 
contribution, dotted line for the MSSM correction with = 0°, dashed hne for the 
MSSM correction with 7(+)7(+) polarization and cj) = 45°, dash-dotted line for the 
MSSM correction with 7(— )7(— ) polarization and 4> = 45°; (b) the cross section 
(t(±^) as a function of sohd line for tree level contribution, dotted line for the 
MSSM correction with (p = 0°, dashed line for the MSSM correction with = 45°; 
(c) the CP-violating parameter ^cp ^ ^ function of V§. 

Fig.3 With m~g = 5 GeV, mi^ = 50 GeV and mi^ = 250 GeV: (a) the relative 
corrections to the cross section of the tt production process a{±±) as a function of 
v^, dotted hne for the MSSM correction with = 0°, dashed hne for the MSSM 
correction with 7(+)7(+) polarization and 4> = 45°, dash-dotted line for the MSSM 
correction with 7(— )7(— ) polarization and cj) — 45°. (b) the relative corrections to 
the cross section cr(±^) as a function of v^, dotted line for the MSSM correction 
with 0=0°, dashed line for the MSSM correction with = 45°; (c) the CP- violating 
parameter as a function of V§. 
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Fig.4 With nif^ = 100 GeV, rrif^ = 450 GeV and = 500 GeV: (a) the cross 
section a of the tt production process from unpolarized photon coUision as a function 
of rrig, sohd hne for tree level contribution, dotted hne for the MSSM correction with 
= 0°, dashed line for the MSSM correction with — 45°; (b) the CP- violating 
parameter ^(^p as a function of rrig. 

Fig.5 With m-g = 100 GeV, = 450 GeV and VI = 500 GeV: (a) the cross 
section a of the tt production process for unpolarized photon collision as a function 
of m^^ , sohd line for tree level contribution, dotted line for the MSSM correction with 
(f) — 0°, dashed line for the MSSM correction with (p — 45°; (b) the CP- violating 
parameter as a function of mf^. 

Fig. 6 With rrig = rrif^ = 100 GeV and \/§ = 500 GeV: (a) the cross section 
(7 of the tt production process from unpolarized photon collision as a function of 

, sohd line for tree level contribution, dotted hne for the MSSM correction with 
= 0°, dashed line for the MSSM correction with — 45°; (b) the CP- violating 
parameter ^cp as a function of mf^. 
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